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ABSTRACT 

Using simulations of slowly rotating stratified turbulence, we show that the a effect responsible for the 
generation of astrophysical magnetic fields is proportional to the logarithmic gradient of kinetic energy density 
rather than that of momentum, as was previously thought. This result is in agreement with a new analytic theory 
developed in this paper for large Reynolds numbers. Thus, the contribution of density stratification is less 
important than that of turbulent velocity. The a effect and other turbulent transport coefficients are determined 
by means of the test-field method. In addition to forced turbulence, we also investigate supernova-driven 
turbulence and stellar convection. In some cases (intermediate rotation rate for forced turbulence, convection 
with intermediate temperature stratification, and supernova-driven turbulence) we find that the contribution of 
density stratification might be even less important than suggested by the analytic theory. 
Subject headings: magnetohydrodynamics (MHD) - Sun: dynamo - turbulence 



1. INTRODUCTION 

Turbulent dynamos occur in many astrophysical sit- 
uations. They tend to develop large-scale magnetic 
structures in space and time that are generally under- 
stood in terms of mean-field dynamo theory (e.g. [M offatt 
19781; lParkerlll979l: IKrause & Radlerl 119801: IZeldovich etafl 
19831 iRuzmaikin et all 119881: IRiidiger & Hollerbachl 12004 
Brandenburg & Subramanian 2005). Central to this theory is 
the a effect, which is a pseudo-scalar multiplying the mean 
magnetic field to give a contribution to the mean electromo- 
tive force. Such a pseudo-scalar can be the result of rotation, 
J7, combined with stratification of density and/or turbulence 
intensity, Vp and/or V« rms , respectively. Here, p is the gas 
density and u rms = (it 2 ) 1 / 2 is the rms value of the turbulent 
velocity, u. 

There have been a number of analytic studies quanti- 
fying the effects of rot ating stratified turbulence on the 
mean electromotive force (|Krau se & Radler 1980; Kichatinc 
19911: IRiidiger & Kichatinovl Il993t iKitchatinov et all 1199 
Radler et al. 2003; Kleeorin & Rogachevskii 120031) . n par- 
ticular, it was found, using the quasi-linear approach (or 
second-order correlation approximation), that the diago- 
nal components of the a tensor for slow rotation rate (or 
small Coriolis number s) are given by (ISteenbeck et al.lfl966h 
IKraus e & Radlel [l980l) 



;f2 • Vln(pw rms ), 



(1) 



where £ a — TQU lms is a relevant length scale and tq is the 
characteristic turbulent time related to the turnover time. 

In the solar convective zone the fluid density p changes by 
seven orders of magnitudes, while the turbulent kinetic en- 
ergy \pu 2 changes by only three orders of mag nitudes and, 
according to mixing length theory ( Vitense 1953), the energy 
flux pu 3 would be approximately constant in the solar con- 
vective zone. This issue has become timely, because there is a 
new numerical technique that allows the different proposals to 
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be examin ed with sufficient accuracy . The so-called test-field 
method ( Schrin ner et alj |2005, 2007) allows one to determine 
all the relevant turbulent transport coefficients in the expres- 
sion for the mean electromotive force without the restrictions 
of some of the analytic approaches such as the quasi-linear ap- 
proach, the path-integral approach or the r approach. With the 
test-field method one solves sets of equations for the small- 
scale fields resulting from different prescribed mean fields — 
the test fields. These equations resemble the usual induction 
equation, except that they contain an additional inhomoge- 
neous term. 

This method is quite powerful because it has been shown 
to be rather accurate and it gives not only the tensor coef- 
ficients of a effect and turbulent diffusivity, but it also al- 
lows the scale-dependence to be determined, which means 
that these coefficients are actually integral kernels that al- 
low the effects of neighboring points in space and time to 
be t aken into account. For det ails r egarding scale separation, 
see iBra ndenbur get alJ ([2008) and Hubbard & Brand enburg! 
(2009). We apply this method to numerical simulations of 
forced turbulence in a stably stratified layer in the presence 
of rotation and a prescribed vertical dependence of the turbu- 
lence intensity. We also use test-field method in simulations 
of turbulent convection and supernova-driven turbulence of 
the interstellar medium (ISM). 

The goal of the present paper is to determine the correct 
scaling of the a effect with density and rms velocity for slow 
rotation, i.e., when the Coriolis number Co = 2Q,tq is much 
less than unity, and large Reynolds numbers. In addition to 
the parameter £ a , we determine the exponent a in the diagonal 
components of the a tensor, 

a = -£l(n- V)\n(p°u lms ). (2) 

Such an ansatz was also made by Riidi ger & Kichatinovl 
(1993), who found that in the high conductivity limit, a = 
3/2 for slow rotation and a = 1 for rapid rotation. However, 
as we will show in this paper, both numerically (for forced 
turbulence and for turbulent convection with stronger tem- 
perature stratification and overshoot layer) as well as analyti- 
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cally, our results for slow rotation and large fluid and magnetic 
Reynolds numbers are consistent with a = 1/2. However, in 
some simulations we also found a = 1/3. 

2. THEORETICAL PREDICTIONS 

We consider the kinematic problem, i.e., we neglect the 
feedback of the magnetic field on the turbulent fluid flow. 
We use a mean field approach whereby velocity, pressure and 
magnetic field are separated into mean and fluctuating parts. 
Unlike in earlier derivations, and to maintain maximum gen- 
erality, we allow the characteristic scales of the fluid density, 
p, the turbulent kinetic energy, \pu 2 , and the variations of 

pu 3 to be different. The strategy of our analytic derivation is 
to determine the £1 dependencies of the second moments for 
the velocity Uj(t, x) Uj(t, x) and for the cross-helicity tensor 

bi(t,x)uj(t,x), where b are fluctuations of magnetic field 
produced by tangling of the large-scale field. To this end we 
use the equations for fluctuations of velocity and magnetic 
field in rotating turbulence, which are obtained by subtracting 
equations for the mean fields from the corresponding equa- 
tions for the instantaneous fields. 

2.1. Governing equations 

The equations for the fluctuations of velocity and magnetic 
fields are given by 



at p 



(3) 



db 



: (B • V)« - (u ■ V)B - B{V ■ u) + Af(b), (4) 



where Equation © is written in the reference frame rotating 
with constant angular velocity ft and we consider an isother- 
mal equation of state, p = c 2 p. The density scale height is 
then constant. Here, c s = const is the sound speed, B is 
the mean magnetic field, p and p are the fluid pressure and 
density, respectively. The terms Af(u) and A/"(6), which in- 
clude nonlinear and molecular viscous and dissipative terms, 
are given by 



Af(u) = (u ■ V)u - (u ■ V)« + f v (u), 



Af(b) = V x [uxb-uxb — 77 V x 6) 



(5) 
(6) 



where pf v (it) is the molecular viscous force and r\ is the mag- 
netic diffusion due to the electrical conductivity of the fluid. 
The fluid velocity satisfies the continuity equation. 

2.2. The derivation procedure 

To study rotating turbulence we perform derivations which 
include the following steps: 

(i) use new variables (V, H) and perform derivations 
which include the fluctuations of rescaled velocity V = p a u 
and the mean magnetic field H = B/p a ; 

(ii) derive equations for the second moments of the veloc- 
ity fluctuations Vi Vj and the cross-helicity tensor bi Vj in k 
space; 

(iii) app ly the spectral closure, e.g., the spectral r approx- 
imation (Pouqu et et al.|[l976t iKleeorin et al.llT990h for large 
fluid and magnetic Reynolds numbers and solve the derived 
second-moment equations in k space; 

(iv) return to physical space to obtain formulae for the 
Reynolds stress and the cross-helicity tensor as functions of 



f2. This allows us to obtain the £1 dependence of the a effect. 
In these derivations we only take into account effects that are 
linear in both A and ft, where A = —V In p. 

To exclude the pressure term from Equation © we take 
twice the curl of the momentum equation written in the new 
variables. The equations which follow from Equations © 
and ©, are given by: 



dt 



V w (\-V)- (v (A) ) V 



= 2 f n ■ v (A) 



V (A) x V) + 2 (ft x V (A)N ) ( V (A) • V 



db 



N{V) 

(fl.vw)v-(v-vw)5 

+(2a - 1) (V (A) ■V)H + Af(b), 



(7) 



(8) 



V + crA. Furthermore, N(V) = 
is the nonlinear term. We consider two 



where V (A) = 

Vx (VxAf(u 

cases: (i) low Mach numbers, where the fluid velocity fluc- 
tuations V satisfy the equation V • V = (1 — cr){V ■ A) in 
the anelastic approximation, and (ii) fluid flow with arbitrary 
Mach numbers, where velocity fluctuations satisfy to standard 
continuity equation. To derive Equation (O we use the identi- 
ties given in Appendix lAl 

2.3. Two-scale approach 

We apply the two-scale approach, e.g., a correlation func- 
tion, 



J J dki dk 2 Vi(ki)Vj(k 2 ) exp{i(ki-x 
+k 2 -y)} = J J dk dK f i:j (k, K) exp(ifc-r + iK R) 

= J dk fij(k, R) exp(ife-r), 

(see, e.g., Roberts & Soward 1975). Hereafter we omit the ar- 
gument t in the correlation functions, fij(k, R) = L(Vi \ Vj), 
where 



L(a; c) 



a(k + \K)c{-k + \K) exp (iK-R) dK, 



and we have introduced the new variables R = \(x + y), 

r = x — y, K = ki+ k 2 , k = h(ki — k 2 ). The variables R 
and K correspond to large scales, while r and k correspond 
to small scales. This implies that we assume that there exists 
a separation of scales, i.e., the maximum scale of turbulent 
motions £0 is much smaller than the characteristic scale Lb 
of inhomogeneity of the mean magnetic field. 

2.4. Equations for the second moments 

Using Equations (Q and ([H) written in k space 
we derive equations for the following correlation func- 
tions: fij(k,K) 



h{tM)V 3 {t,k 2 ). 
tions are given by 

dfii 



= Vi{t,ki)Vj{t,k-z) and gij(k,K) = 
The equations for these correlation func- 



dt 



u A (1) f ■ 



(2). 



^(V, k f 
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ilqfcqfnj H~ ^npq^pl fiq ^ifq m 



+4+ ATf, 



0t 



-2n 



y ■ 



(9) 



~t~ npg &pZJ 



+ 



iK ( , 



"J 



(10) 



where fcj 



kikj/k 2 , AW = ki — icrA, and similarly for 



A(2), = A^ , Aj i 1 V(A( 1 )) 2 , and similarly for The 
source terms and in Equations (0) and ( TTOb contain 
large-scale spatial derivatives of H, which describe the con- 
tributions to turbulent diffusion (if. and 1^ are gi ven by Equa- 
tions (A4) and (A6) in lRogachevskii & Kleeorinl (120041) ). The 
terms Tv/y and Afgij are related to the third-order moments 
that are due to the nonlinear terms and are given by 



A%=P iro (A(i))Aqy ro (fci)]VS-(fc 2 ) 



-V l (fc 1 )P jm (A( 2 ))AA[T/ m (fc 2 )] ! 



Af 9ij =Af[b l (k 1 )]V ] (k 2 ) + 6 i (fc 1 )P im (A( 2 ))AT[y m (fe2)] ! 
where Py (A) = <Jy — Ay . 

2.5. T-approach 

The equations for the second-order moments contain 
highe r order moments and a cl o sure problem ari ses (lOrszagl 
Il970t iMonin & Yagloml 119751: iMcCombl 119901) . We ap- 
ply the spectral r approximation o r the third-order clo- 
sure procedure dPouquet et"aT] Il976t iKleeorin et al.l 119901 : 
Rogachevskii & Kleeorin 2004). The spectral r approxima- 
tion postulates that the deviations of the third-order-moment 
terms, A//y(fc), from the contributions to these terms af- 
forded by the background turbulence, Mf^\k), are ex- 
pressed through similar deviations of the second moments, 



jy(fc)-/W(fc),i.e 



r r (k) 



(11) 



and similarly for the tensor Afgij . Here the superscript (0) 
corresponds to the background turbulence (i.e., non-rotating 
turbulence with a zero mean magnetic field), r r (k) is the char- 
acteristic relaxation time of the statistical moments, which can 
be identified with the correlation time r(k) of the turbulent 
velocity field for large Reynolds numbers. We also take into 

account that (fe) = 0. We apply the r-approximation ( fTTT i 
only to study the deviations from the background turbulence. 
The statistical properties of the background turbulence are 
assumed to be known (see below). A justification for the 
r approximation in different situations has been obtained 
thro ugh numerical simulations and ana l ytical studies (see, 
e.g., Brandenburg & Subramanian 2005; Rogachevskii et al. 
120111) . 



2.6. Model for the background compressible turbulence 
We use the following model for the background turbulence: 



/, 



^ (0) (W 



(0) 



(fca) = { 



A«A„ (2) 



5 {j AW . A^ 



+i (Xikj - Xjh) + p c ki kj + - (kiS7 



jVlj J/87rfc4(i + Mc / 2 ) 



V 2 + 0(£ 2 X 2 ), (12) 



where 



fi c = (V • it) 2 / (V x it) 2 < 1 



(13) 



is the degree of compressibility of the turbulent velocity field, 
the terms oc p c take into account finite Mach numbers com- 
pressibility effects, r(fc) = 2rof(fc), E(k) = —d?(k)/dk, 
f(k) = (fc/fco) 1-9 , 1 < q < 3 is the exponent of the 
spectrum function (q = 5/3 for Kolmogorov spectrum), 
fco = (-q 1 , io is the maximum scale of turbulent motions, 
and uq is the characteristic turbulent velocity at scale Iq. The 
motion in the background turbulence is assumed to be non- 
helical. For low Mach numbers (/i c <C 1), Equation dT2l > 
satisfies to the condition: V • V = (1 — <r)(V ■ A). 

2.7. Contributions to the a effect caused by rotation 
Since our goal is to determine the a effect, we solve Equa- 
tions (O and ( [Tol l neglecting the sources l(- and If - with large- 
scale sratial derivatives of H . We subtract from Equations (0 
and (flOb the corresponding equations written for the back- 
ground turbulence, and use the spectral r approximation. We 
take into account the effects which are linear in A and in fi. 
We also assume that the characteristic time of variation of the 
second moments is substantially larger than the correlation 
time r(k) for all turbulence scales. This allows us to get a 
stationary solution for Equations (O and ( TTOb for the second- 
order moments. Using this solution we determine the con- 
tributions to the mean electromotive force caused by rotating 
turbulence: 



bi(k)Vj(-k)dk =e 



mji 



£mji 
- H i\ n kpl 



dk 



3i 



ij(k) dk 



p(0) 
/ nq 



(0) 
npqjqj 



( s jpqfl 

x(eM$ +*»p q k jpl fW)]. (14) 
After performing the integration in k space, we get: 

(fiiVj + Qj-Vi - ASijSl ■ V) In V 2 
+ (2ct - 1) (QiXj + CljXi - 46ijSl ■ A) 
+ \p c (QjXi - AiliXj + Sijfl ■ A) . (15) 

For low Mach numbers (p c -C 1) and for a = 1/2, the a 
tensor depends only on V In V 2 , i.e., 



- U 2 - 

15 3 



< 
15 



(fiiVj + fij-Vi - 4%0 • V) In V 2 , (16) 



where we used = ayPj, and ay = §(ay + a^) is the 
symmetric part of ay. Furthermore, the pumping velocity, 
7i = jEi nm a mn is independent of rotation for small Coriolis 
numbers, because the rotational contribution to the pumping 
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velocity vanishes. In this case the pumping velocity is deter- 
mined only by the inhomogeneity of turbulent magnetic dif- 
fusivity. 

For arbitrary Mach numbers and when 



1 . 



(17) 



the diagonal part of the a effect depends only on V In V 2 , 



i.e., 



1(5/2 

15 1 ' 



(18) 



n=-^° (n-V)ln(p CT Wrms ), 



The latter equation can be rewritten in the following form: 

.- » - - ■■ ; i j. i / ( 1 

15 

where cr is determined by Equation ( fTTI i. In the next section 
we will determine the exponent a from numerical simulations. 

3. NUMERICAL SIMULATIONS 

We now aim to explore the universality of the derived scal- 
ing by comparison with results from very different astrophys- 
ical environments. To do this, we performed three quite dif- 
ferent types of simulations of: 

(i) artificially forced turbulence of a rotating stratified gas, 
where the density scale height is constant and the turbulence 
is driven by plane wave forcing with a given wavenumber. 

(ii) supernova-driven interstellar turbulence in a vertically- 
stratified local Cartesian model employing the shearing sheet 
approximation; 

(iii) turbulent convection with and without overshoot layers. 
The advantage of the first method is that it allows to impose 

a well-defined vertical gradient of turbulent intensity, i.e., of 
the rms velocity of the turbulence such that V In u rms is ap- 
proximately constant over a certain z interval, excluding the 
region near the vertical boundaries. Naturally, the physically 
more complex scenarios (ii) and (iii) are less-well controlled 
but allow to demonstrate the existence of the a effect scaling 
in applications of direct interest to the astrophysical commu- 
nity. We use different variants of the test-field method to mea- 
sure the a effect (and other turbulent transport coefficients). 

3.1. Basic equations 

In the following we consider simulations where the turbu- 
lence is either driven by a forcing function /, in which case 
we assume an isothermal gas with constant sound speed c s , or 
it is driven through heating and cooling either by supernovae 
in the interstellar medium or through convection with heating 
from below. In all cases, we solve the equations for the veloc- 
ity, U, and the density p in the reference frame rotating with 
constant angular velocity fl and linear shear rate S: 



p— = -Vp + p(f n + f + g) + V-(2ispS), (20) 



DU 

Dt 
dp 

dt 



- V • P U, 



(21) 



where /q = (2D,U y , (S — 2Sl)U x , 0) is a combined Coriolis 
and tidal acceleration for a linear shear flow Us = (0, Sx, 0), 
v is the kinematic viscosity, / is a forcing function, g is grav- 
ity, and Sij — \ {V jUi + VjC/,*) — • U is the traceless 
rate-of-strain tensor, not to be confused with the shear rate S, 
and D/Di = d/dt+ (U + Us) • V is the advective derivative 



with respect to the total (including shearing) velocity. In the 
isothermal case, the pressure is given by p — pc 2 , while in all 
other cases we also solve an energy equation, for example in 
terms of the specific entropy s = c v In p — c p In p, where c p 
and c v are respectively the specific heats at constant pressure 
and constant volume, their ratio 7 = c p /c v is chosen to be 
5/3, and s obeys 



PT 



Ds 
Dt 



-VF rad + /9 r-p 2 A + 2^ j oS 2 , (22) 



where the temperature T obeys (c p — c v )T = p/p, F la a is 
the radiative flux, and T is a heating function. In the isother- 
mal case, the entropy equation is not used, and the forcing 
function / consists of random, white-in-time, plane, non- 
polarized waves with a certain average wavenumber, kf . 

The simulations are performed with the PENCIL Code 
(http://pencil-code.googlecode.com) which uses sixth-order 
explicit finite difference s in space and a third-order a ccurate 
time stepping method (Brandenburg & Doblerl 12002). The 
simulations of supernova-driven turbulence in the ISM have 
been performed using the NlRVANA-III code (Ziegleri r2004l) 
with explicit viscosity and resistivity. 

3.2. The test-field method 

We apply the kinematic test-field method (see , e.g., 
iSchrinneret ail 12005. 2001: iBrandenburg et al.l l2008) to de- 
termine all relevant turbulent transport coefficients in the gen- 
eral relation 

Si = otijBj + r]ijkBj t h, (23) 

where Bj.k = VfeBj is the magnetic gradient tensor. The 
test-field method works with a set of test fields B T , where 
the superscript T stands for the different test fields. The cor- 
responding mean electromotive forces £ are then £ T , and are 
calculated from £ T — u x b T , where b T = V x a T with 



0a T 



I h 11 f! 11 b 1 )' + r/VW , (24) 



72 _T 



dt 



with U and u taken from the solutions of the momentum 
equation. In the case with shear, we replace da T /dt by 
da T j dt + Us ■ Va T + Sa^x. On the top and bottom bound- 
aries we assume perfect conductors, and for the x and y direc- 
tions periodic boundary conditions. These small-scale fields 
are then used to determine the electromotive force £ T corre- 
sponding to the test field B T . The number and form of the 
test fields used depends on the problem at hand. 

We either use planar (xy) averages, which depend only on z 
and t (hereafter referred to as test-field method I), or, alterna- 
tively, we assume that the mean field varies also in the x and 
y directions, but that the turbulence is homogeneous in those 
two directions, and that the z direction constitutes a preferred 
direction of the turbulence (test-field method II). In the former 
case only the x and y components of £ are important for dy- 
namo action, and the magnetic gradient tensor has only two 
non-vanishing components which can be expressed in terms 
of the components of the mean current density, J = V x B, 
alone. Here and elsewhere, we adopt units for which the mag- 
netic permeability is unity. Thus, we have 



£% Q^ij-^j Vij J j 7 



(25) 



with i and j being either 1 or 2. 
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Alternatively, the mean electromotive force is assumed to 
be characterized by only one preferred direction which we 
describe by the unit vector e. Then, £ can be represented 
in the form 

£ = aj_B + (an — a±)(e ■ B)e + je x B 

—r]±J — (r]\\ —r)±)(e ■ J)e — Se x J (26) 
—kj_K — (km — K±)(e ■ K)e — pe x K 

with nine coefficients a±, a\\, . . ., p. Like J = V x B, 
also K is determined by the gradient tensor VB. While J is 
given by its antisymmetric part, K is a vector defined by K = 
e-(VB) s with (Vi?) s being the symmetric part of VB. For 
details of this method, referred to below as test-field method 
II, see |Brandenburg et al.l (120121) . Equation (O is expected to 
apply to a±, while ai i is know to behave differently and can 
have opposite sign ( Brandenburg et al. 1990; Ferrierel ll992t 
Riidi ger & Kichatinovlll993l) . 

Errors are estimated by dividing the time series into three 
equally long parts and computing time averages for each of 
them. The largest departure from the time average computed 
over the entire time series represents an estimate of the error. 

3.3. Simulations of forced turbulence 
We begin by studying forced turbulence. We consider 



a domain of size L 7 



Li, 



L z in Cartesian coordinates 



(x,y,z), with periodic boundary conditions in the x and y 
directions and stress-free, perfectly conducting boundaries at 
top and bottom, z = ±L z /2. The gravitational acceleration, 
g = (0, 0, — g), is chosen such that the density scale height 
H p = Cg I g is small compared with the vertical extent of the 
domain L z . The smallest wavenumber that fits into the cu- 
bic domain of size L 3 is k\ — 2ir/L, so the density contrast 
between bottom and top is exp(27r) w 535 and the density 
varies like p = poexp(—z/H p ), where po is a constant. In 
all cases, we use a scale separation ratio kf jk\ = 5, a fluid 
Reynolds number Re = u rms /vkf between 60 and 100, and a 
magnetic Prandtl number Pr^/ = v jr\ of unity. In our units, 
c s = k\ = po = /j,o = 1. We use a numerical resolution of 
128 3 mesh points for all forced turbulence runs. 

We perform simulations for different values of the rms ve- 
locity gradient, \' u ) = dlnu rms /dz, but fixed logarithmic 
density gradient, dlnp/dz = — A. Thus, we have 



e 2 a nx (<j-a w /a) 



(27) 



i.e., a can be obtained conveniently as the value of A*- 11 -* / A 
for which a vanishes. By arranging the turbulence such that 
A("> and A are approximately independent of z, the value of a 
is also approximately constant. In that case, however, all the 
other turbulent transport coefficients are z-dependent. How- 
ever, by normalizing 7 by u rms /6 and the other coefficients 
by r]to{z) — ii rms (z)/3fc/, we obtain non-dimensional quan- 
tities that are approximately independent of z. We denote 
the corresponding non-dimensional quantities by a tilde and 
quote in the following their average values over an interval 
z\ < z < Z2, in which these ratios are approximately con- 
stant. 

In FigureQ] we plot the normalized profiles of a± and a>\\ as 
functions of z. Note that there is a limited range where both 
functions are approximately constant. For A*") ^ 0.5, they 
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FIG. 1. — Comparison of a± (solid) and an (dashed) for the runs with 
A(") = 0.33 and 0.61. 

are of opposite sign; see the lower panel of Figure Q] and Ta- 
ble Q] As discussed above, this be havior has been seen and 
interpreted in earlier calculations dBrandenburgetalJll990l 
lFerrierdll992l:lRudiger & Kichatinovll 19931) ! 
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\( u ) f or Co = 0.15 (upper panel) and 0.4 (lower panel), giving respectively 
cr « 1 /2 and 1 /3 as the zeros in each graph. 

In Figure|2]we plot the dependence of the normalized mean 
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values of a± on the value of for two values of Co. The 
value of a can then be read off as the zero of that graph. We 
find cr » 1/2 for low values of Co, and a somewhat smaller 
value (a « 1/3) for larger values of Co. 
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FIG. 4. — Results for the other transport coefficients Co = 0.15. 



FIG. 3. — Transport coefficients for the run with A("> = 0.61. 

In Figure [3] we show the z-dependence of the remaining 7 
normalized coefficients. They are all approximately indepen- 
dent of z within the range z\ < z < z-x. In Figure [4] we 
show the scaling of these turbulent transport coefficients (in- 
cluding now also au) with \( u \ The results for an and 7 
suggest a dependence proportional to the gradient of p ai u lms 
with <j a . between and 0.3 for a\\ and <r 7 between 0.7 and 
1 for 7. On the other hand, all the other coefficients seem to 
be independent of A and we find /3\\ « /3± w 1.1, 8 w 0, 
K11 w —0.1, k± rs —0.4, and /2 « —0.05. These results 
are quantitatively and qualitatively in agreement with those 
of iBrandenb urg et al] (12012b . The fact t hat S turned out to 
be es sentially zero was addressed earlier ( IBrandenburg et aO 
2012), where it was found that significant values are only 
found for scale separation ratios around unity, i.e., when the 
scale of the mean field is comparable to that of the turbulent 
eddies. 

3.4. The case of supernova-driven ISM turbulence 

We now tu rn to simulations of supernova-driven ISM tur- 
bulence (see iGressel et al.1 l2008allbL for a detail e d des crip- 
tion of the model ), similar to those of iKorpi et al.l ( 1 19991) and 
iGent et aT] (12012b . but extending to larger box sizes, thus al- 
lowing for better scale separation. In these simulations, ex- 
pansion waves are driven via localized injection of thermal 






z [kpc] 



t=1 ,90Myr, ct = 0.327 
red. x 2=1 - 01 

2 4 



FIG. 5. — Time-averaged vertic al profile of from ISM turbulence (see 
Gressel. Elstner & Rudiger 201 1), obtained with the TF method (gray line). 
The best-fit model (black line) with a = 0.327 is obtained by method of 
least squares, weighted with the standard deviation (shaded area) in ot^ . 



energy, Tg^a;, t). Additionally, optically thin radiative cool- 
ing with a realistic cooling function A = A(T) and heating 
pT{z) lead to a segregation of the system into multiple ISM 
phases. 

Here we solve the visco-resistive compressible MHD equa- 
tions (supplemented by a total en ergy equation), using the 
NlRVANA-III code dZieglerl 120041) : for the full set of equa- 
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TABLE 1 

Basic parameters and turbulent transport coefficients for the forced turbulence simulations. 



Ke A/ 


K, U 


A v * 


rv , 1-2 /Ok. 

Of 1 liftj 












K±/ r lto( z ) 


K \\/Vto( z ) 




70 


0.14 


0.13 


0.47 ±0.16 


0.1 ±0.1 


-0.1 ± 0.1 


1.5 ±0.2 


2.1 ±0.3 


-0.2 ±0.0 


-0.1 ± 0.1 


-0.6 ±0.1 


-1.8 ±0.4 


56 


0.17 


0.25 


0.82 ±0.10 


0.7 ±0.1 


0.0 ±0.1 


1.1 ±0.1 


1.0 ±0.1 


0.0 ±0.0 


-0.2 ± 0.0 


-0.7 ± 0.1 


-0.1 ±0.1 


62 


0.16 


0.33 


0.57 ±0.03 


0.5 ±0.1 


-0.1 ± 0.1 


1.1 ±0.2 


1.2 ±0.2 


0.0 ±0.0 


-0.1 ± 0.0 


-0.6 ±0.1 


-0.3 ±0.2 


47 


0.20 


0.46 


0.32 ±0.26 


0.2 ±0.3 


-0.2 ± 0.1 


1.1 ±0.2 


1.0 ±0.2 


-0.0 ±0.0 


-0.1 ± 0.0 


-0.4 ± 0.0 


-0.1 ±0.0 


65 


0.15 


0.61 


-0.42 ±0.06 


0.1 ±0.5 


-0.4 ±0.0 


1.1 ±0.1 


1.2 ±0.1 


-0.1 ±0.1 


-0.1 ± 0.1 


-0.4 ± 0.1 


-0.3 ±0.2 


59 


0.16 


0.97 


-1.29 ±0.26 


-0.2 ±0.5 


-0.6 ±0.2 


0.8 ±0.3 


0.7 ±0.3 


-0.1 ±0.0 


-0.1 ± 0.1 


-0.3 ±0.1 


0.1 ±0.1 


108 


0.36 


0.10 


0.20 ±0.02 


-0.1 ±0.1 


0.0 ±0.0 


0.9 ±0.1 


1.6 ±0.2 


0.1 ±0.0 


0.1 ± 0.0 


-0.5 ±0.1 


-1.6 ±0.2 


93 


0.41 


0.14 


0.12 ±0.04 


-0.1 ±0.0 


-0.0 ±0.0 


1.0 ±0.1 


2.2 ±0.4 


0.0 ±0.0 


0.1 ± 0.0 


-0.3 ±0.1 


-1.7±0.4 


75 


0.51 


0.35 


0.30 ±0.05 


0.3 ±0.0 


-0.2 ± 0.1 


1.0 ±0.1 


0.8 ±0.1 


0.0 ±0.0 


-0.3 ±0.0 


-1.1 ± 0.0 


0.2 ±0.0 


74 


0.52 


0.59 


-0.36 ±0.02 


0.2 ±0.1 


-0.7 ±0.1 


1.2 ±0.0 


0.8 ±0.0 


-0.3 ±0.1 


-0.2 ± 0.1 


-0.9 ±0.0 


0.3 ±0.0 
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FIG. 6. — Likelihood map based on the reduced-x 2 error estimate, with the 
best-fit parameter set indicated by the cross. 



tions, we refer the reader to Equation (2.1) in lGresseH (2010). 
For the present run, we chose a resolution of 128 x 128 x 512 
mesh points, and apply a value of Pr^/ = 2.5. The fluid 
Reynolds number, defined as Re = u rms £q/{2ttv), varies 
within the domain and takes values Re ~ 70-165, which is 
somewhat larger than for the forced turbulence case. The 
same applies to the Coriolis number Co = 2f2ro which is 
here ~ 0.24. 

Despite the mentioned differences, the basic properties of 
the turbulence producing an a effect are in fact quite similar: 
rotation of the system together with stratification in the mean 
density and turbulence amplitude. Notably, the ISM simula- 
tions are strongly compressible with peak Mach numbers of 
up to ten, corresponding to a typical value of /i c ~ 1.9, and 
with peak values up to five. We note that the ISM simula- 
tions also include shear with S = — 57, which may affect the 
production of vorticity. 

To obtain an estimate for a from the time-series of a sin- 
gle simulation run, we apply a method distinct from the one 
described above: We treat the a ( p ( p(z) profile (inferred with 
test-field method I) as the data to be modeled and obtain er- 
ror estimates by means of the standard deviation within four 
equal sub-intervals in time (see gray line and shaded areas 
in Figure [5). We then compute time averages of the profiles 
for V In p and V In u rms (here without error estimates), from 
which we compute a model prediction for a ( j >( j > (z) based on 
the expression 

«<W> = -^"L fi ■ V In (p a u rms ) , (28) 



where we have assumed that l\ can be replaced by 

£ a — 7" a u rms (z) , (29) 

with u lms (z) being the z-dependent rms velocity and r a (as- 
sumed independent of z) is a characteristic timescale related 
to the a effect. We apply a least-square optimization allow- 
ing r a and a as free parameters. The best-fit model according 
to Equation d28l is plotted as a black line in Figure [5] along 
with the stated values for r a and a, and matches well the data 
within the error bars. 

Because we cannot a priori assure that r Q is uniform in 
space, and because this might affect the precise determination 
of (T, we perform an additional test. We do this by indepen- 
dently estimating £o(z) from the two-point velocity correla- 
tion function, computed in horizontal slabs around a given 
galactic height z, and time averaged over multiple snapshots. 
By comparing the obtained £q(z) with u lms (z), we find that 
our data are broadly compatible with a uniform To of about 
1.2 Myr. To corroborate the fit, we compute a likelihood map 
in the parameter space spanned by r Q and a and find that the 
best-fit parameter set is, in fact, located at the global mini- 
mum of the reduced-x 2 map; see Figure [6] The best-fit value 
of r a is around 1.9 Myr, which is indeed compatible with the 
correlation time tq. 

To conclude this section, we remark that we here find a 
somewhat smaller exponent, a ~ 1/3, which suggests that 
this case deviates from the theoretical prediction. It has how- 
ever a similar exponent as in our stratified forced turbulence 
simulations with larger values of Co. Note that with the 
determined value for /i c , Equation (F7\ predicts a value of 
a ~ 0.62, which is a factor two larger than obtained from 
the fit. The reason for such discrepancy between the theoret- 
ical predictions and the simulations of supernova-driven ISM 
turbulence might be caused by the fact that the theory is devel- 
oped for simplified conditions which are different from these 
simulations. 

3.5. Convection-driven turbulence 

Many astrophysical bodies have turbulent convection 
zones. Again, rotation and stratification induce helicity into 
the flow and therefore drive an a effect. Using mixing length 
theory (I Vitensd 1 1 95 3h, we know that the con nective flux is 
proportional to pu^ ms (Brandenbur g et alJ2005l) . In the steady 
state, the total energy flux is constant in space, so if most of 
the flux is carried by convection, then pu^ ms — const and 
thus its vertical gradient vanishes. If the scaling of Section [3T4l 
were applicable also to this case, i.e., if a cx dlnpw,^ ms /dz, 
then a would vanish. 
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FIG. 7. — Dependence of a on z for convective turbulence without over- 
shoot layers (gray line, with shaded areas indicating fluctuations) compared 
with Equation {2} (black line / dotted) applying a = 0.76 and l a = 0.41 d, 
as obtained from a least-square fit within the highlighted interval in z/d. 
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FIG. 8. — Same as Figure[7j but the case with overshoot layers and weak 
temperature stratification, comparing a (gray line, shaded areas) measured 
via the TF method with Equation (2) (black line / dotted) yielding a = 0.35 
and i a = 0.6 d as best-fit values. 
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FIG. 9. — Same as Figure|7] but the case with overshoot layers and stronger 
temperature stratification, comparing a (gray line, shaded areas) with Equa- 
tion {2) using a = 0.48 and i a = 0.56 d (black line / dotted). 

To investigate this somewhat worrisome possibility, we now 
consider a simulation of turbulent convection in a stratified 



TABLE 2 

Summary of convection runs. 



Run 



Ra 



Ta 



Ap, 



total 



Apcz overshoot Res. 



A 1.3-10° 
B 6.1 ■ 10 5 
C 6.1 ■ 10 5 



3.6 ■ 10 J 
6.4 ■ 10 4 
6.4 ■ 10 4 



64 
37 
290 



64 
7 
296 



256 J 
128 3 



128 J 



layer, heated from below by a constant energy flux z ■ -F X ad 
at the bottom, where we adopt the diffusion approximation 
for an optically thick gas with -F ra d = — A'VT and radia- 
tive conductivity K. We either consider a constant value of 
X = K/ (pc p ) with an enhanced turbulent heat conductivity 
Yt nea r the surface as in spherical simulations of Kapy la et alJ 
d201 ll) . or, alternatively, a piecewise constant profile K(z), 
such that g ■ Vs is positive in the middle of the domain, which 
corresponds to co nvective instability. The latter setup is de- 
scribed in detail in Kapyla et al.l (120091) . The hydrostatic equi- 
librium value of g- Vs is proportional to the Ray leigh n umber, 
Ra, which is here around 10 6 ; see lKapyla et alJ J2008) for the 
definition. The rotational influence is here measured in terms 
of the Taylor number Ta = (2Cld 2 /v) 2 , which is 3.6 • 10 5 
for Run A. Summary of our convection simulations is given 
in Table |2] Run A without overshoot layers and 6.4 • 10 3 for 
Runs B and C. Here d is the thickness of the unstable layer. 
The density contrast in the run without overshoot layers is 64. 
The vertical boundary conditions are stress-free and we use 
Pr/ij = 0.5 and Re^/ = 52 in Run A, and Pi~m = 1 and 
ReM = 13 in Runs B and C. In Figure [7] we give the results 
for Run A without overshoot layers. The Coriolis number 
varies with height and is about 0.2 in the middle of the layer. 
The rotation axis is anti-parallel to the direction of gravity, 
corresponding thus to a location at the north pole. The sys- 
tem is therefore isotropic in the xy plane and we consequently 
quote the mean between the two horizontal components, i.e., 
a = h(a xx + ctyy), using test-field method I, which corre- 
sponds to a± of test-field method II. Note that in this sim- 
ulation, a(z) shows a sinusoidal profile, suggestive of either 
weak stratification or effects of boundaries. Given that den- 
sity stratification is not small (factor 64), it is plausible that 
the effects of boundaries are here responsible for the extended 
regime with negative a. 

Our simulation shows that the best fit is obtained for a ps 
0.75. A possible reason for this unexpected behavior might 
be poorer scale separation in convection simulations com- 
pared with forced turbulence simulations. The other possi- 
bility is related to the absence of convective overshoot layers 
discussed above. This idea is partly confirmed by compar- 
ing with simulations that include convective overshoot layers. 
Now the best fit value for a is found to be about 1/3 when 
the temperature stratification is weak (Figure[8]l and about 1/2 
when the temperature stratification is strong (Figure[9]). 

4. CONCLUSIONS 

While the present investigations confirm the old result 
that the a effect in mean-field dynamo theory emerges as 
the combined action of rotation and stratification of either 
density or of turbulent intensity, they also now point to- 

TABLE 3 
Summary of results for <j. 





Co 


Ap/p 


CT 


forced turbulence 


0.15 


535 


1/2 




0.40 


535 


1/3 


supernova-driven ISM 


0.24 


1000 


1/3 


convective turbulence (CT) 


0.2 


64 


3/4 


CT with overshoot 


0.2 


37 


1/3 




0.2 


290 


1/2 


analytic theory 






1/2 
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ward a rev ision of the standard f ormula for a. The old 
formula by ISteenbeck et al.l {T966) predicted that the effect 
of stratification can be subsumed into a dependence on the 
gradient of pu r ms- This for mula was then generalized by 
Rudiger & Kichatinov (1993) to a dependence on p CT it rms , 
where a — 3/2 in the high conductivity limit for slow ro- 
tation, and (7 = 1 for faster rotation. In contrast, our new 
results now clearly favor a value of a below unity. The ide- 
alized case of artificially forced turbulence can most directly 
be compared to our analytic derivation, since it agrees in all 
the made assumptions. The obtained value of u = 1/2 agrees 
very well with the theoretical expectation. A similar expo- 
nent is found for the case of turbulent convection with higher 
temperature stratification, but the results seem to depend sen- 
sitively on model parameters (see Table [3]). Here more de- 
tailed studies will be required. Moreover, the result a = 1/2 
arises naturally from analytical considerations for large fluid 
and magnetic Reynolds numbers and slow rotation as the only 
tenable choice, but those considerations have not yet been per- 
formed for the cases of intermediate and rapid rotation. 

Forced turbulence simulations show a trend toward smaller 
values of a around 1/3 for faster rotation and also in cases 
of supernova-driven turbulence. Turbulent convection with 
overshoot also gives 1/3 in one case of moderate temperature 
stratification with overshoot, while simulations without over- 
shoot point toward values somewhat larger values around 3/4. 
However, in none of the cases have we found that the a ef- 
fect diminishes to zero as a result of a trend toward constant 
convective flux where puf ms is approximately constant. 



In spite of the considerable scatter of the values of a found 
from various simulations, it is worth emphasizing that in all 
cases a is well below unity. On theoretical grounds, the value 
1/2 is to be expected. Except for the forced turbulence simu- 
lations that also yield 1/2 for slow rotation, all other cases are 
too complex to expect agreement with our theory that ignores, 
for example, inhomogeneoities of the density scale height and 
finite scale separation. 
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APPENDIX 

A. IDENTITIES USED FOR THE DERIVATION OF EQUATION Q 

To derive Equation (O we use the following identities: 

2" 



v! A) vf-%(v( A )) 



Vj, 



P° [Vx(Vx4 = 

p a [vx[vx(uxn)]] 
x(vxn)j, 

v W> - s* (v (A) ) 2 ] (vxnjj = (nxvW^A-v) 
+{n.v w ){v (x) xv) i . 

Equation ( IA3b is obtained by multiplying the identity 

by A 2 V}, where is the fully antisymmetric Levi-Civita tensor, A mn = A m A n /A 2 , and the identity ( IA4b is valid for arbitrary 
vectors f2 and A. 



(Al) 

(A2) 

(A3) 
(A4) 
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